Cold Fermionic Atoms in Two-Dimensional Traps — Pairing versus Hund's Rule 
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The microscopic properties of few interacting cold fermionic atoms confined in a two-dimensional 
(2D) harmonic trap are studied by numerical diagonalization. For repulsive interactions, a strong 
shell structure dominates, with Hund's rule acting at its extreme for the mid-shell configurations. 
In the attractive case, odd-even oscillations due to pairing occur simultaneously with deformations 
in the internal structure of the ground states, as seen from pair correlation functions. 

PACS numbers: 03.75.Hh, 05.30.Fk, 45.50.Jf, 67.85.-d 



Quantum dots and related nanostructures appeared as 
the first "artificial atoms" [1], where semiconductor tech- 
niques allowed to design the quantum confinement of a 
few particles in the laboratory (see [2] for a review). In 
more recent years, artificial quantum confinement also 
became possible with ultracold trapped atomic gases Q : 
Magneto-optical techniques allow to create very clean 
systems where the trapped particles can be bosons as 
well as fermions. Strength and sign of their interac- 
tions can be tuned, making them repulsive or attractive 
via Feshbach resonances, allowing for a variety of quan- 
tum many-body phenomena to be studied theoretically 
as well as experimentally. Cold atoms may be loaded into 
"deep" optical lattices, where the single sites may resem- 
ble harmonic confinement. With cold, trapped fermions, 
it became possible to experimentally study the crossover 
from a weakly attractive Fermi gas to the unitary limit 
of infinite scattering length, and beyond Q. Attractive 
fermionic atoms in traps share many of their properties 
with nuclei, such as Bardeen- Cooper- Schrieffer (BCS) 
pairing and the occurence of shell structure [Hj]. 

The "Coulomb blockade" in the discrete charging and 
discharging of quantum dots leads to significant conduc- 
tance oscillations. These are manifest in the "fundamen- 
tal energy gap" A 2 (TV) = A X (TV + 1) - Ai(TV), given 
by the difference between the modulus of the particle- 
removal and the addition energy, i.e., by the difference 
in the chemical potential of a quantum system confining 
TV + 1 and TV particles, Ai(iV) = E (N) -E (N-1), 
where Eq(N) is the TV-body ground state energy. Very 
recently, Cheinet et al. [6] reported on the realization of 
a similar interaction blockade [7| in the transport of cold 
atoms through a double well in an optical lattice, observ- 
ing discrete steps in the well population with increasing 
"bias" potentials [note that the atom number distribu- 
tion as a function of the bias relates to A2(TV)]. These 
findings indeed bring implementations of "atomtronics" , 
as suggested in the pioneering work of Seaman et al. [8[ , 
much closer to realization. 

In this Letter, we thus investigate effects of interaction 
blockade for quantum dot-like, 2D fermionic atom traps. 
For repulsive interactions, a strong shell structure occurs 



and Hund's rule acts at its extreme at half-filled shells. 
For attractive interactions, however, a pronounced odd- 
even staggering in the ground state energies (as recently 
discussed for 3D systems by Zinner et al. [9]) is accom- 
panied by broken symmetries in the internal structure of 
the many-body states. 

Let us now consider TV atoms of spin 1/2 and equal 
masses m confined in a 2D harmonic trap, interacting 
through a contact potential, 
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where g' is the coupling constant. We use as energy unit 
hwQ and as length unit i = (h/mujo) 1 / 2 . The dimension- 
less coupling constant is g = g f /(huo£ 2 ). The few-body 
problem of Eq. (pQ) is solved by means of the full con- 
figuration interaction (CI) method [1 Qj ] - In more than 
one spatial dimension, contact interactions do not allow 
to diagonalize the many-body Hamiltonian unless regu- 
larized properly. An approximate way to achieve this is 
a momentum-cutoff to a given subspac e effectively 
renormalizing the interaction strength [12[. In the ap- 
plied model space of the six lowest oscillator shells we 
can then relate the interaction strength g to the 2D scat- 
tering length a by comparing the CI two-body ground- 
state ener gy w ith the exact result derived by Busch and 
coworkers [13| , as shown in Ref. LL2|. The obtained values 
of a are listed in note fill . 

Expanding the interacting TV-body state on the given 
set of Slater determinants and diagonalizing [151 ] . the CI 
method provides the many-body energies and wave func- 
tions of both ground states and (low-lying) excitations, 
which are eigenstates of the total spin S and z-projection 
of the orbital angular momentum M . Figure [T] shows the 
fundamental gap A2(TV) vs TV. In the upper panel, for 
repulsive interactions g > 0, one clearly recognizes the 
distinctive features of a shell structure^ strikingly similar 
to those observed for quantum dots [16[ (confirming the 
earlier prediction by Capelle et al. [3])- The two peaks for 
TV = 2, 6 form when the first two energy shells of the 2D 
harmonic oscillator are completely filled by atoms as they 
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FIG. 1: (Colour) A 2 (A) vs N. (a) The green, red, blue, black 
curves are calculated by using g = 0.3, 1, 3, 5. (b) The green, 
red, blue, black curves are calculated by using g = -0.3, -1, 
-3, -5. (c) Results from the seniority model. 



are consecutively loaded into the trap. Since the height 
of the peaks is mainly dictated by the energy spacing 
(huoo = 1) between neighboring shells, A2(N = 2,6) w 1. 
A2(N) depends only weakly on the interaction strength 
g, as it increases from g = 0.3 (green curve) up to g = 5 
(black curve). The peak for N = 4 at mid-shell between 
N = 2 and N = 6 is a consequence of Hund's rule acting 
at its extreme, where up to half-filling of the degenerate 
shell the Pauli principle eliminates the interactions. 

The signatures of attractive interaction (g < 0) are 
dramatically different from the repulsive case of Fig.QJa). 
Figure [2(b) shows that the non- interacting shell structure 
is gradually washed out as the magnitude of atom-atom 
attractive interaction is increased. As g goes from g = 
—0.3 (green curve) up to g = —5 (black curve), A2(N) 
develops a pronounced even-odd pattern, with A2 > 
(A2 < 0) for even (odd) A", corresponding to S = 
(S — 1/2). Pairing effects are likely to be the dominant 
mechanism behind these odd-even oscillations. If N is 
even, all atoms form singlet pairs and a large amount 
of energy is required to add one unpaired atom. For 
odd N, A2(N) < since energy is gained by pairing 
with an opposite-spin particle. However, the odd-even 
staggering of Fig. [2(b) could originate from either pairing 
or Jahn- Teller effects, as discussed for nuclei, metallic 
grains, clusters [171], and even for 2D systems [18[. 
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FIG. 2: (Colour) E int /N 3 / 2 vs N for g = ±0.3. 



In order to isolate effects from pairing we study the 
seniority model from nuclear physics (also called the 
method of 'exact pairing' [19]), applied to the 2D sys- 
tem. The model restricts the interactions in the 2D trap 
to be exclusively of the pairing type, that act only be- 
tween particles in time-reversed orbits. Unpaired par- 
ticles do not participate in the interaction beyond the 
mean-field level, except for blocking certain final states 
available for pairing. Using an attractive ^-function pair- 
ing potential as for the full microscopic model above, we 
truncate the model space after the fourth oscillator shell. 
The results are shown in Fig. [2(c), comparing well with 
the full many-body calculations. For weak interactions, 
g = —0.3, A 2 (A^) displays the non-interacting shell struc- 
ture. For stronger g, a pronounced odd-even staggering 
appears. We find large gaps between the ground and 
first excited state in even-, and the absence of a gap 
in odd- numbered systems. To discern between inter- and 
intra-shell pairing [5] , we study the fraction of the pairing 
energy that comes from the off-diagonal Hamiltonian ma- 
trix elements, which involve pairs being moved between 
the oscillator shells. For smaller interaction strengths, 
this energy is a fraction of a percent. Increasing the cou- 
pling to g = —5, the off-diagonal contribution is larger 
than 20%. This was roughly the same for all particle 
numbers, though partially filled shells have higher con- 
tributions than particle numbers around closed shells. 

Figure [2] shows the interaction energy E[ nt for 2 < 
N < 20 for small g [20[, where the non-interacting con- 
tribution was subtracted from the total energy, E[ nt = 
Eq(N) — Eo(N,g = 0). In order to magnify the fine 
structure of E{ nt in Fig. [2] we replace it with the scaled 
quantity Ei nt /N 3 / 2 , since in the Thomas- Fermi approxi- 
mation Ei nt ~ A/ -3 / 2 . The plot clearly shows that in the 
case of repulsive interaction (g = 0.3) the energy gain is 
minimum for closed shells (N = 2, 6, 12, 20). Conversely, 
the energy gain is largest for the half-filled open shells at 
Af = 4, 9, 16. In fact, in those cases the spin is maximized 
due to Hund's rule (S = 1, 3/2, 2, respectively). In the 
attractive case, already for weak interactions the pattern 
of E[ nt /N 3 / 2 shows a marked even-odd alternation: the 
energetically favored ground states occur at even N when 
all atoms are paired. 

The key quantity for g < is the pairing energy gap 
A, measuring the interaction-energy gain (expense) by 



M = M = 2 M = 4 




FIG. 3: (Colour) Excitation energies vs M for N = 8. The 
top, middle, bottom panels are obtained for \g\ = 0.3, 3, 5, 
respectively. The red (blue) circles label the case of attractive 
(repulsive) interaction. Only the six lowest excited states are 
shown for each value of M. The filled circles mark the ground 
states, whose energies are taken as reference. 



adding an atom to the TV-body system if N is odd (even). 
We define A as A = [Ai(iV) - Ai(7V + l)]/2 with N 
odd. This is plotted in Fig. [4] vs g for the open-shell 
cases N = 3, 7. Remarkably, the black (N = 7) and red 
(N = 3) curves almost overlap, showing that pairing is a 
generic feature of the few-body system. 

Complementary information on the pairing energy A 
is obtained by the analysis of the excitation spectrum for 
fixed even N. Figure [3] focuses on the low-lying energy 
levels of N = 8 as the magnitude of g increases, going 
from \g\ = 0.3 (top panel) up to \g\ = 5 (bottom panel). 
We resolve the six lowest-energy levels for orbital angular 
momenta < M < 4. For small values of g (g = ±0.3, 
top panel) the non-interacting low-lying sequence clearly 
emerges: three M = states, two M = 2 states, and one 
M = 4 state may be obtained by variously arranging two 
atoms with opposite spin in the third open shell of the 2D 
harmonic oscillator, while the remaining six atoms fill the 
two lowest shells. To have the lowest M = 1 or M = 3 
excitations one particle (or hole) must be excited into 
the 4th (2nd) shell, whereas to have higher-energy states 
with even M, one particle (hole) must be excited into 
the 5th (1st) shell or two into the 4th (2nd) shell. As 
\g\ increases (center and bottom panels of Fig. [3]), two 
fundamentally different types of excitation spectra ap- 
pear for positive (blue curves) and negative (red curves) 
g. For g < 0, all excited states appear at some large en- 
ergy related to the pairing gap A: A pair must be broken 
to excite the system, no matter the value of M. As \g\ 
increases, the gap becomes uniform with respect to M, 
indicating the onset of intershell pairing. The opposite 
holds for g > 0, since excitation energies decrease for in- 
creasing g, implying that several excited states become 




FIG. 4: (Colour) Alternative estimates of the pairing energy 
gap A vs g for N = 4 (3), 8 (7). 




FIG. 5: (Colour) Plot of Pn( r i r o) in the xy plane. The black 
dot at (x,y) = (xo, 0) locates the spin- J, reference atom (xo is 
the average value of r). Top row: g = 3. Bottom row: g — —3. 
From left to right column, N = 2, 4, 6, 8, respectively. The 
squares' size is 5 x 5, and 15 equally spaced contour levels go 
from blue (minimum) to red (maximum). 

almost degenerate with respect to the ground state, as in 
strongly correlated quantum dots [2]. 

The lowest excitation energy, E exc (N), may be used 
as an alternative definition of A [21]. Such quantity is 
plotted in Fig. E] for N = 4 (blue line) and N = 8 (green 
line) and may be compared with the previous evaluations 
of A through the chemical potential. We see that all the 
different estimates basically coincide, independent of N. 
A fair agreement is found also for cases other than open 
shells, provided \g\ is sufficiently large. It is interesting to 
compare the curves for E exc (N) in Fig. 2] with first-order 
perturbation theory in g (magenta curve). At small \g\ all 
estimates of A agree with the predicted value of — g/(47r), 
which is linear in g. Around g ~ — 1, however, strong 
deviations from the perturbation- theory result arise, due 
to the non-linear behavior of A. 

We finally investigate the internal structure of ground 
state wave functions, computing the conditional proba- 
bility P crcro (r, ro) of finding an atom with spin a = j, j at 
r provided another one is fixed at ro with spin g$: 

Pa* (r,r ) = A^2(6(r - ri)S a(Tz S(ro - rj)5 aoCTj ) , (2) 

where A is a normalization factor. The probability 
P|l(r, ro) is plotted in Fig. [S] for even TV, with ro = (xq, 0) 
labeled by a black circle [for convenience, we scaled 
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the maximum height of P[i(r, Vq) to the same value in 
each panel] . Figure [5] only shows correlations between 
opposite-spin atoms, since the parallel spin probability 
P aocro (r, ro) is almost structureless, except the exchange 
hole around the fixed atom (the contact interaction does 
not scatter atoms with parallel spin). 

Already for N = 2 (left column of Fig. [5]), the difference 
between the repulsive (top panel, g = 3) and attractive 
(bottom panel, g = —3) case is manifest. In fact, whereas 
for g > the |-atom may be found in the antipodal po- 
sition with respect to the fixed j-atom, for g < both 
atoms tend to overlap in space, suggesting the formation 
of a bound pair. The contour plots for g = — 3 (bottom 
row) show a peak of the probability of finding an j-atom 
close to the j-atom, and they are rotated by 2tt/N with 
respect to the repulsive case, as displayed in the upper 
panel. Intriguingly, whereas the overall spatial distribu- 
tion for P|j,( r ' r o) tends to an isotropic distribution of 
atoms for g = 3, for g = — 3 the distribution is strongly 
distorted in space. This distortion is particularly clear 
for N = 4 and N = 6, with six particles filling a shell 
in the non-interacting limit, which should maintain cir- 
cular symmetry. The origin of the deformation seen in 
Pjj might be attributed to the Jahn- Teller effect or the 
arrangement of bound pairs in the trap, as discussed by 
Stecher et al [22[ for the four-particle system. 

In summary, we studied shell structure and energy gaps 
for a few fermionic particles with repulsive as well as 
attractive contact interactions in a 2D harmonic trap. 
For repulsive interactions, shell structure and Hund's 
rule lead to a significantly enhanced fundamental energy 
gap at closed-shell and mid-shell configurations. For at- 
tractive interactions on the BCS-side of unitarity, how- 
ever, a pronounced odd-even staggering was found, where 
the ab-initio results agree well with the seniority model. 
Interaction blockade, as discussed here for 2D "quan- 
tum dots" with cold atoms, may experimentally be ob- 
served in atom-transport studies, as recently performed 
by Cheinet and coworkers [6[ for bosons, and an exten- 
sion to the fermionic case is called for. 
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